Stable Quantum Resonances in Atom Optics 
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A theory for stabilization of quantum resonances by a mechanism similar to one leading to 
classical resonances in nonlinear systems is presented. It explains recent surprising experimental 
results, obtained for cold Cesium atoms when driven in the presence of gravity, and leads to further 
predictions. The theory makes use of invariance properties of the system, that are similar to those 
of solids, allowing for separation into independent kicked rotor problems. The analysis relies on a 
fictitious classical limit where the small parameter is not Planck's constant, but rather the detuning 
from the frequency that is resonant in absence of gravity. 

PACS numbers: 05.45-Mt, 03.75.-b, 42.50.Vk 



(N 

Q 

U 
d 



> 

o 
o 
in 
o 

<N 
O 



> 

x 



Resonances are among the most typical manifestations 
of the quantum mechanical behavior. They are charac- 
terized by a precise relation between the parameters of 
the system and an external driving; for instance, reso- 
nant excitation of an atom requires that the difference 
between two levels be equal to the frequency of the driv- 
ing field times Planck's constant. Classical resonances 
of nonlinear systems are more robust. Regions of ap- 
preciable size in phase space may respond in a resonant 
way, while other regions are chaotic. In recent experi- 
ments it was found that a quantum system that is 
far from its semiclassical limit exhibits resonant behavior 
very similar in nature to the one found for classical non- 
linear systems: it is robust, and the high precision that 
is characteristic of the usual quantum resonances is not 
required. It is observed for various values of the param- 
eters and of the initial conditions. In this Letter the ex- 
perimental results are explained and further predictions 
are made. Details are given in [Q. The theory reveals 
how resonances, similar to those characteristic of classi- 
cal systems, may arise in a quantum system, far from the 
semiclassical limit. The theoretical model of the experi- 
ment is a variant of the kicked rotor, that is a standard 
system used in the investigation of classical Hamiltonian 
chaos and its manifestations in quantum mechanical sys- 
tems In experimental realizations of the model, 
laser-cooled atoms are driven by application of a stand- 
ing electromagnetic wave. The frequency of the wave is 
slightly detuned from resonance, so a dipole moment is 
induced in the atom. This moment couples with the driv- 
ing field, giving rise to a net force on the center of mass of 
the atoms, proportional to the square of the electric field 
HI). As the wave is periodic in space, the atom is thus 
subjected to a periodic potential. The wave is turned 
on and off periodically in time, and the time it is on is 
much shorter then the time it is off. A realization of a 
periodically kicked particle is then obtained. For atoms 
driven in the horizontal direction, effects of gravity are 



negligible, and localization in momentum is observed [Bj, 
as theoretically expected (see below). 

In the experiments that provide the subject of the the- 
oretical analysis of the present letter, laser cooled Cesium 
atoms were driven in the vertical direction, and gravity 
was found to produce new remarkable effects. Localiza- 
tion was destroyed; furthermore, in the vicinity of the 
frequencies that are resonant in absence of gravity, a new 
type of ballistic spread in momentum was experimentally 
observed |t|-[| . A fraction of the atoms were steadily ac- 
celerated, at a rate which was faster or slower than the 
gravitational acceleration depending on what side of the 
resonance the driving frequency was. Such atoms were 
exempt from the diffusive spread that took place for the 
other atoms, and their acceleration depended on the dif- 
ference between the driving and the resonant frequencies. 
The dynamics of the atoms that are falling as a result of 
gravity and are kicked by the external field is modeled 
by the time-dependent Hamiltonian: 
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where t' is the continuous time variable, the integer vari- 
able t counts the kicks, P, X are the momentum and the 
position operator respectively. Units are chosen so that 
the mass of the atoms is 1, the Planck's constant is 1, and 
the spatial period of the kicks is 2tt. The dimensionless 
parameters k,r),r fully characterize the dynamics. They 
are expressed in terms of the physical parameters as fol- 
lows: k = k/H, t = hTG 2 /M, T) = MgT/(hG), where 
M, T, k, g are the mass of the atoms, the kicking period, 
the kick strength, and the gravitational acceleration re- 
spectively, and 2ir/G is the spatial period of the kicks. 
The positive a;— direction is that of the gravitational ac- 
celeration. Throughout the following, time is a discrete 
variable, given by the kick counter t. 

It is expedient to measure the momentum in the free 
falling frame, notably to replace P — -i' by P, resulting 
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in the time-dependent Hamiltonian : 
H(t') = hp+Ht') 2 + kcos(X) S(t'-tr). (2) 
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FIG. 1. Phase portraits for the map (W) on the 2-torus, 
with n/r = 0.01579. Left: stable (1,0) fixed point J = 0, 
0o = 0.42, for t = 5.86, k = 1.329. Right: stable (10,1) peri- 
odic orbit at r = 6.31, k = 0.067. 

This Hamiltonian is related to ([!]) by the gauge trans- 
formation e ir > xt '/ T , The main advantage of (||) compared 
to (Q) is that it depends on X only via cos(X). Conse- 
quently the evolution only mixes momenta which differ 
by integers: hence, quasi-momentum is conserved. In 
absence of gravity (77 = 0), (g) reduces to the Hamilto- 
nian of the kicked rotor (classically resulting in the Stan- 
dard Map) with stochasticity parameter K = kr. For 
K > K c « 0.9716 diffusion in momentum takes place. 
For values of K near to integer multiples of 2ir accelerator 
modes are found Q . These are stable phase space islands 
that travel ballistically, resulting in linear, rather than 
diffusive growth of momentum with time. The quantum 
mechanical study of the systems described by (^) starts 
with decomposing the momentum as p — n + /3, where 
n, (3 denote the integer and the fractional parts of the 
momentum p respectively. Since the quasimomentum j3 
is conserved only n varies in the course of the dynamics. 
The evolution is equivalent to that of a superposition of 
independent kicked rotors, each characterized by a dif- 
ferent value of (3. Such a rotor will be called /3-rotor and 
the one step evolution operator is: 
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where 9 = x mod(27r) and the momentum operator is 
N = —i-ja- For (3 — and 77 = this is the usual 
kicked rotor. The classical diffusion is then suppressed 
by a mechanism that is similar to Anderson localization 
in disordered solids P,|To|]. For r = 2nl/m, where / and 
m are integers, the eigenstates of (|J) are extended in mo- 
mentum and ballistic growth of momentum takes place 
in most cases. This is the quantum resonance |Tl| ]. The 
Talbot length corresponds to 47r. For fixed (3 ^ and for 
typical irrational t/2tt the classical diffusion is again sup- 
pressed, as was observed in experiments on laser-cooled 



Sodium and Cesium atoms || . Quantum resonances only 
occur at special values of f3, hence their experimental ob- 
servation is difficult [^2) , since it is impossible to prepare 
all atoms with the same quasimomentum (3. Averaging 
over (3 results in linear (rather than quadratic) growth of 
the squared momentum [Q. This demonstrates the high 
sensitivity of the quantum resonances to fine experimen- 
tal details. 

In the presence of gravity, 77 ^ 0, the evolution operator 
(^) is time dependent. The localization in momentum is 
destroyed, as it was experimentally observed jl|-[| . While 
this is expected for time dependent random potentials, 
the present time dependence is just quasi-periodic, so a 
deeper analysis is required. However, the most surprising 
experimental result was that an appreciable fraction of 
atoms were found to accelerate (in the free falling frame) 
for various values of the experimental parameters, for 
values of r in intervals of appreciable size around integer 
multiples of 27r. Here is a quantum resonance that is ro- 
bust, in contrast to the usual quantum resonances. It is 
reminiscent of classical nonlinear resonances, namely the 
accelerator modes of the Standard Map; however, it has 
no counterpart in the classical limit of (^j. In what fol- 
lows a theoretical explanation of this effect is presented 
§• 

We consider the case when r is close to a resonant 
value 2irl (I > integer), and the kicking strength k is 
large. We hence write r = 27r/ + e, k = k/\e\ with e small. 
Noting e~ llxln = e ~ i7ri ", (||) takes the form (apart from 
an irrelevant phase factor): 



Up® 



jkcosO ^— -^j H p (I ,t) 



where 



\e\N: 



Hp{I, t) = \ sign(e)/ 2 + I(nl + t{(3 + tr) + g)) 
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If |e| is assigned the role of Planck's constant, then (|^J^) 
is the formal quantization of either of the following clas- 
sical (time-dependent) maps : 



h+i = h + fcsin( 



8 t+1 = 9 t ±I t + nl + T (J3 +trj + rj/2) mod(27r), 1 ' 

where ± has to be chosen according to the sign of e. The 
small |e| asymptotics of the quantum [3— rotor is thus 
equivalent to a quasi-classical approximation based on 
the "classical" dynamics @. We emphasize that "clas- 
sical" here is not related to the h — > limit but to the 
limit e — > instead. The two limits are actually incom- 
patible with each other except possibly when I = 0. For 
the sake of clarity the term "e— classical" will be used in 
the following. 

Changing variable to J t = It ± irl ± t(/3 + rjt + rj/2) 
removes the explicit time dependence of the maps (|(j), 
yielding: 
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Jt+i = Jt + ~k sin(0 t+ i) ± TTj , O t +i =6 t ±J t . (7) 

If J, 6* are taken mod(27r), then (^) define maps of the 
2— torus in itself. Let Jo, 8q be a period p-fixed point of 
either of the toral maps thus defined. Then iteration of 
(0) yields, at t = p : 



J P = Jo + 27rj , 9 p = 9 Q + 2irm 



(8) 



for some integers j , to. In terms of the original dynamics 
@ this yields a family of orbits such that, for all integer 
t, 

P t = 6»o = ^o mod(27r) , I tp =I + apt , (9) 

where a = T T '7 + 27rj /p and 

I = J TnlTT((3 + r)/2) + 2nm' , (10) 

with to' any integer. Thus primitive periodic orbits of the 
toral maps (Q) correspond to families of accelerator orbits 
of the dynamics @, marked by linear average growth of 
momentum with time. If the periodic orbits are stable, 
then the accelerator orbits are surrounded by islands of 
positive measure in phase space, also leading to ballistic 
(linear) average growth of momentum in time. These are 
named accelerator modes. They are characterized by the 
integer couple (p,j) formed by the order p and by the 
jumping index j. 
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FIG. 2. Contour plots of Husimi functions at times 

t =0,2,4,8,16 for a f3— rotor initially prepared in the coher- 
ent state centered at the e— classical (1,0) accelerator mode 
of Fig.l (left). The black spots in the centers of the contours 
are an ensemble of classical points of size h = 1, evolving 
according to the e— classical dynamics (Q). 

Modes of order 1 correspond to fixed points of (Q). 
They are given on the 2-torus by Jo = 0, #o = where 



sin(0j) = (27Tj =F rr])/k 



(11) 



and j is any integer such that the rhs is not larger than 1 
in absolute value. A linear stability analysis shows that 



for any integer j each map (fjj) has exactly one stable 
period- 1 fixed point on the 2-torus, given by ( |ll"| ) if, and 
only if, 



|27Tj =p TTJ | < k < y/16 + (27TJ =F TT]) 2 



(12) 



When k increases beyond the right-hand limit, such fixed 
points turn unstable and bifurcations occur. At small 
values of k higher-order accelerator modes appear, asso- 
ciated with higher-period primary orbits. Examples of 
stable periodic orbits near r = 27r are presented in Fig.l. 
Initial physical momenta no = |e| _1 /o for e— classical ac- 



celerator modes are obtained from (10) for any < /3 < 1. 
If the stable islands associated with e— classical acceler- 
ator modes have a large area compared to |e|, then they 
may trap some of the rotor's wave packet and give rise to 
quantum accelerator modes traveling in physical momen- 
tum space with speed ~ a/\e\ = — rrj/e + 27rj/(p|e|). In 
order that such modes may be observed, the phase space 
distribution associated with the initial rotor state must 
significantly overlap the islands. This picture is con- 
firmed by numerical simulations. Fig. 2 shows the quan- 
tum phase-space evolution of a (3— rotor with j3 = 0.2188 
started in the coherent state centered at the position of 
a (1,0)- accelerator mode. The Husimi functions com- 
puted at subsequent times closely follow the motion of 
the e— classical mode. In the e— semiclassical regime, 
accelerator modes are expected to decay exponentially 
in time due to quantum tunneling out of the classical is- 
lands, with decay rate 7 e oc exp — (const. /|e|). This decay 
was also numerically confirmed [Q. 

In the experiments described in , the initial state 
of the falling atoms is satisfactorily described by an inco- 
herent Gaussian mixture of momentum eigenstates cen- 
tered at p = 0, with rms deviation ~ 2.55 in our units. 
This is equivalent to a statistical ensemble of /3-rotors. 
The above theory describes this situation for small e, 
where the e- classical approximation holds. From (||) it 
is found that the e-classical (pj)-mode started at t = 
with Iq = no|e| is located at time t at the momentum: 



n — tiq — tr-q/e + 27rfj/(p|e| 



(13) 



In Fig. 3 this prediction is compared to the results of nu- 
merical simulations in the vicinity of r = 27r for k = 0.87T 
and t]/t — 0.01579 (that is the value of the gravitational 
acceleration in our units). The time is t = 60 kicks. The 
initial state is an incoherent mixture of 50 momentum 
states, distributed as reported in experimental papers. 
Lines are as predicted by (|l3|). Good agreement is found. 
The mode (1,0) was identified experimentally (as can be 
seen on comparing Fig. 3 to Fig. 2 of 0; also recall that 
our positive direction is that of the gravitational acceler- 
ation). Further higher order modes are predicted in the 
present work. Longer time is required for their resolu- 
tion, because they move slower, yet they may be eventu- 
ally resolved, see the inset of Fig. 3. Similar behavior is 
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predicted for r w 47r (as described in Q), and for higher 
multiples of 2~k as well. We finally remark that (i) quan- 
tization of the /3— rotors momentum enhances quantum 
modes at discrete values of (3, given by @ with Iq an 
integer multiple of e. This results in preferred values 
for physical momentum p = n + f3, spaced by « 2irl/r. 
This was experimentally detected, and it is seen in the 
fine structure of the modes in Fig. 3; (ii) the accelerator 
modes of the /3-rotors result in acceleration of the par- 
ticle erent 
supe 
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FIG. 3. Momentum distribution in the falling frame at 
t = 60 kicks for different values of the kicking period r near 
t — 2ir. Note the negative sign of p. Darker regions cor- 
respond to higher probability. The initial state is a mixture 
of plane waves sampled from a gaussian distribution of mo- 
menta. Full lines are the theoretical curves ([[3]), with orders 
and jumping indexes as indicated by the arrows. The inset 
shows data at t = 400. 



The key steps in the theory presented in this work 
are: (i) Transformation to the free falling frame leading 
to the Hamiltonian (») Conservation of the quasi- 
momentum 0, allowing for introduction of the /3-rotors; 
(Hi) Approximation of the motion by the one generated 
by the e-classical Hamiltonian (j^) near t = 2irl; (iv) e- 
classical accelerator modes (|l0|) dominate the dynamics 
near r = 2irl. 

What stabilizes these accelerator modes? On top of 
the long wavelength behavior related to the resonance 
t = 2ttI, there is a short wavelength behavior. Its dy- 
namics is dominated by the fictitious classical mechanics 
termed e-classical. This is reminiscent of the "caustics 
without rays" [^3) . The generality of this mechanism for 
nonlinear stabilization of quantum resonances should be 
explored. The fact that the intermediate-time dynamics 



is dominated by a discrete set of modes, which exponen- 
tially decay in time, bears a distinct resemblance to the 
Wannier-Stark problem of a Bloch particle in a constant 
field p^ j. How far this analogy carries is, in our opin- 
ion, an interesting theoretical question. Other theoretical 
questions well worth further analysis are about the na- 
ture of the asymptotic (in time) behavior (in particular, 
the mechanism for destruction of localization), and the 
nature of the e-semiclassical systematic expansion. On 
the experimental side, it would be interesting to classify 
the various accelerator modes, namely to perform "ac- 
celerator mode spectroscopy" ; and also to experimentally 
explore differences in the quantum dynamics between co- 
herent and incoherent superpositions in the initial state 
of the atoms. This may be of importance for atomic in- 
terferometry and for quantum computation. 
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